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Abstract. Very high energy gamma-rays from blazars traversing cosmological distances through the metagalactic 
radiation field can convert into electron-positron pairs in photon-photon collisions. The converted gamma-rays 
initiate electromagnetic cascades driven by inverse-Compton scattering off the microwave background photons. 
Using a model for the time-dependent metagalactic radiation field consistent with all currently available far- 
infrared-to-optical data, we calculate the cascade contribution from faint, unresolved high- and low-peaked blazars 
to the extragalactic gamma-ray background as measured by EGRET. For low-peaked blazars, we adopt a spectral 
index consistent with the mean spectral index of EGRET detected blazars, and the EGRET luminosity function. 
For high-peaked blazars, we adopt template spectra matching prototype sources observed with air-Cherenkov 
telescopes up to 30 TeV, and a luminosity function based on X-ray measurements. The low number of ~ 20 for 
nearby high-peaked blazars with a flux exceeding 10 -11 cm~ 2 s _1 above 300 GeV inferred from the luminosity 
function is consistent with the results from air-Cherenkov telescope observations. Including the cascade emission 
from higher redshifts, the total high-peaked blazar contribution to the observed gamma-ray background at GeV 
energies can account for ~ 30%. 
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1. Introduction 

Isotropic, diffuse background radiation presumably due 
to faint, unresolved extragalactic sources has been ob- 
served in nearly all energy bands. The confirmation 
of an extragalactic gamma-ray background by EGRET 
(Energetic Gamma-Ray Experiment Telescope) on board 
the Compton Gamma Ray Observatory has extended the 
spectrum up to an energy of ~50 GeV. A first analysis 
of the data resulted in a total flux of (1.45 ± 0.05) • 10 -5 
photons cm~ 2 s _1 sr _1 above 100 MeV and a spectrum 
that could be fitted by a power law with an spectral in- 
dex of -2.10 ± 0.03 (Sreekumar et al 1998). These val- 
ues are strongly dependent on the foreground emission 
model which is subtracted from the observed intensity to 
obtain the extragalactic residual. Since, using the fore- 
ground model by Hunter et al (1997), a residual GeV halo 
remained after subtraction (in addition to the isotropic 
extragalactic background), the foreground model had to 
be improved. This led to a new analysis of the EGRET 
data, and a new result for the extragalactic background 
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spectrum, now showing a dip at GeV energies and an over- 
all weaker intensity of (1.14 ± 0.12) • 10~ 5 photons cm~ 2 
s _1 sr _1 (Strong, Moskalenko & Reimer 2004). This new 
result can help us to understand the origin of the extra- 
galactic background radiation. 

Since EGRET detected a large number of extragalactic 
gamma-ray sources belonging to the blazar class of AGN, 
a reasonable assumption is that the gamma background is 
produced by unresolved AGN. Stecker & Salamon (1996) 
were able to explain 100% of the background by includ- 
ing the effect of variabiliy. But the resulting number of 
faint, nearby blazars was too low compared to the numbers 
detected by EGRET (Chiang & Mukherjee 1998). Using 
a gamma-ray luminosity function from EGRET data, 
Chiang & Mukherjee (1998) found that only 25% to 50% of 
the gamma-ray background could be explained by blazars. 
This result was questioned by Stecker (2001) who argued 
that a statistically independent analysis is contrary to his 
assumption of a correlated radio and gamma-ray emission 
and introduces a bias. Giommi et al. (2006) estimate the 
gamma-ray background contribution of blazars from deep 
blazar counts in the radio band, from multi-frequency sur- 
veys, and from interpolations of data using synchrotron- 
sclf-Compton models. They found that blazars possibly 
can contribute 100% in the cx 0.5 — 50 MeV band. Since 
their calculation over-predicts the observations at ener- 
gies (E > 100 MeV) by a large factor they conclude that 
the duty cycle of blazars must be rather low, and this is 
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consistent with the results from long-term monitoring o 
TeV blazars such as Mrk421 or Mrk501. Due to the gap ir 
sensitivity in the 100 keV to 100 MeV band, the spectra 
shape of blazars between medium energy gamma rays anc 
high energy gamma rays is poorly constrained by obser 
vations. 

Here extending the existing models we assume a pop- 
ulation of BL Lacertae objects (BL Lac) with a spec- 
tral energy distribution such that their flux at EGRET 
energies is generally too low to be detected, while thei] 
very high energy gamma ray flux is strong. Most of these 
sources are at redshifts high enough for pair attenuatioi 
with the UV-IR metagalactic radiation field (MRF) tc 
take place. A significant part of their VHE emission is re 
processed by inverse Compton cascades (Protheroe 1986 
Maciolek-Niedzwiecki, Zdzarski & Coppi 1995, Protheroi 
& Stanev 1993). For a low extragalactic magnetic fielc 
strength (< 10 -17 G) the emission contributes to the 
source spectrum (Dai et al. 2002, Fan, Dai & Wei 2004). 
Assuming a much higher magnetic field of 10 -9 G the 
electrons and positrons become isotropized. The photons 
produce the so-called gamma-ray halo (Aharonian, Coppi 
& Volk 1994). The inverse Compton emission of blazars 
can contribute to the GeV gamma-ray background (Coppi 
k Aharonian 1997). 

The paper is organized as follows: in the next sec- 
tion we present the set of equations used to calculate the 
gamma-ray background including the component from un- 
detected sources and, assuming no extragalactic magnetic 
field, from cascade emission. In Sections 3 and 4 we give 
the choice of parameters for two populations of AGN, the 
undetected EGRET blazars and high peaked BL Lacs. 
The results are shown in Section 5, followed by a discus- 
sion in Section 6. 

In this paper we use a Hubble constant of Hq = 71km 
s -1 Mpc -1 and a flat Universe with the cosmological pa- 
rameters J7m=0.3 and Oa=0.7. 

2. Gamma-ray background 

The emission in units of [GeV -1 cm -2 s -1 sr -1 ] of a popu- 
lation of unresolved gamma-ray sources to the extragalac- 
tic gamma-ray background can be described by 
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with dN 1 /dE^(z s , L) as the intrinsic gamma-ray flux and 
the cascade emission dN c / dE 1 (z s , L) from a source at red- 
shift z 8 , tacitly assuming that the cascade contributes at 
the redshift of the source which is the case for all pair cre- 
ation optical depths t 77 3> 1. This assumption is valid for 
most of the absorbed spectra, which we assume to gener- 
ically extend to 30 TeV, since the majority of sources at 
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Fig. 1. Spectral indices of blazars in the 3rd EGRET cat- 
alog (mean index 1.37 ± 0.04) 
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1 have 7Vv ^> 1 above 50 GeV. The cosmo- 



logical volume element is dV c /dz, L m is the total lumi- 
nosity of the weakest source, dN/ (dV dL) the luminosity 
function and Nq is the normalization constant for a single 
source spectrum defined as a function of the total intrinsic 
luminosity L 
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Here the limits E m i n and E max depend on the definition of 
the energy range of the luminosity used in the luminosity 
function. 



2.1. Extragalactic gamma-ray absorption 

The high energy photons (> 20 GeV) of cosmic gamma- 
ray sources undergo interactions with the UV-FIR meta- 
galactic radiation field (MRF) . Therefore gamma-ray pho- 
tons are absorbed due to pair-production (Kneiske et al. 
2004 and references therein). 

We use a broken power-law for the intrinsic spectrum 
in the energy range from 100 MeV to 30 TeV 
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with ot = oti for EZ, > E' peak and a = a 2 for E^ < E' peak . 

The two slopes and peak energy will be determined by 
the spectrum of the observed sources. 



2.2. Cascade contribution 

The following analytical approximation, similar to the one 
in Fan, Dai & Wei (2004), used to calculate the particle 
flux of a pair-cascade following the pair production process 
from a gamma-ray source at redshift z s . The cascade term 
in Eq. §\§ accounts for the first generation of electrons pro- 
duced in photon-photon pair production. Corrections for 
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the second generation of electrons is negligible since the 
contribution is very small for the chosen gamma energy 
range. For all spectra we assume a maximum energy of 
30 TeV, so the inverse Compton scattering takes place in 
the Thomson regime only. The choice of this energy limit 
is motivated by the highest energies observed in BL Lac 
spectra. Note that the effect of a higher maximum energy 
has been studied in Protheroe & Stanev (1993), and can 
generally be assumed to be weak for plausible spectral in- 
dices a > 2 above 30 TeV. We ignore magnetic fields in 
the following calculations, although in a realistic scenario 
the cascade emission will likely be isotropizcd due to the 
deflections of the electrons and positrons by intergalactic 
magnetic fields. Taking magnetic fields into account, we 
would have to extend the luminosity function to include 
the unbeamed host sources of blazars, the radio galaxies. 
Their greater number density would compensate for the 
lower luminosity resulting from the absence of relativistic 
beaming. We can thus approximate the cascade compo- 
nent given by 
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with f(x) = 2xlog(a;) + x + 1 - 2x 2 , (0 < x < 1) and 
x = (E' 7 / A^e) . The isotropic photon density of the cos- 
mic microwave background v(e) is a black body with a 
temperature of T = (2.728 ± 0.008) (1 + z s ) 
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In this calculation positrons will be treated as elec- 
trons. On average, each absorbed gamma-ray photon pro- 
duces two electrons with an energy e e = 1/2E} / , since 
pair creation has no preferred direction in the center- 
of-mass frame. (Note that we use the index 1 to show 
that the gamma energy is different from the gamma en- 
ergy in Eq. and ([5]).) The electron spectrum dN e /dj e 
as a function of the Lorenz factor j e = e e /mc 2 for 
Ye.min < 7c < 7c, max, 7c, max = 30 TeV/(2m e c 2 ) and 
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=max[i( J Bi/e e ) 1 / 2 , 100 MeV/(2m e c 2 )] is given by 
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For numerical reasons the time integral has been re- 
placed by an inverse-Compton time-scale tic multiplica- 
tion. A comparison between the result using Eq. (J4j)- (J7j) 
and a more exact calculation as in Protheroe & Stanev 
(1993) shows good agreement. 

3. Undetected EGRET blazars 

The luminosity function of resolved EGRET sources, ex- 
tended to the faint end, has been computed by Chiang & 
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Fig. 2. Template spectra for a source of case 1 (thin lines) 
and case 2 (thick lines) located at different redshifts. For 
both cases the intrinsic spectrum (dashed lines), the ab- 
sorbed spectrum (solid lines) and the cascade emission 
(dot-dashed lines) is shown. In the lower right panel the 
data are shown for Mkn501 (see Ref. of Fig.l and Kneiske 
et al. 2004). The shaded region is the intrinsic region for 
different MRF models (Kneiske et al. 2004). In the lower 
left panel the same is shown for 1ES1959+650 at a redshift 
of 0.47. 

Mukherjee (1998). We used their model changing only the 
spectral index from a = 2.10 to a = 2.37. The new spec- 
tral index was determined by fitting the newly determined 
EGRET background spectrum at energies E 1 < 1 GeV. 
This choice is supported by averaging the spectral indices 
for detected blazars from the 3rd EGRET catalog (see 
Figure[T]) which led to a mean spectral index of 2.37±0.04. 
We also included the cascade emission of the sources. 
Although the total intensity of the newly determined spec- 
trum of the extragalactic background has become lower, 
the total contribution of this blazar population (consisting 
mainly of flat-spectrum radio quasars) to the extragalactic 
background in the energy range 100 MeV< E 1 <20 GeV 
amounts to about 60%. 



4. High peaked BL Lacs 

The remaining excess of the measured gamma-ray back- 
ground could be produced by the cascade emission of high- 
energy peaked or X-ray blazars belonging to the HBL 
(XBL) class (def. for HBL see e.g. Ghisellini et al 1998). 
To test this idea we adopt an average HBL spectrum and a 
TeV-luminosity function obtained from X-ray properties. 



4.1. Template spectra 

A number of extragalactic gamma-ray sources have been 
detected with imaging air-Cherenkov telescopes (e.g., Wei 
2006). Some of them were bright enough to resolve their 
spectra in the TeV energy band. At first glance the 
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Table 1. Parameters of HBL template spectra. 
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Fig. 3. Same as in Figure [2] but the flux is shown for a 
source of case 3. The data and intrinsic regions are for 
Mkn421, PKS2005-489 and PKS2155-304. (see refs Fig.l 
and Kneiske et al. 2004). 

spectral energy distributions show a variety of shapes. 
However, by calculating their intrinsic spectra taking the 
extragalactic absorption into account, the spectra in the 
source frames are consistent with a remarkably similar 
shape which can well be approximated with a broken 
power-law. The spectral index above the peak is close to 
~ 2.5 (cf. Mannheim 1998), while it is much flatter below 
the peak (< 1.7 in order to remain below the EGRET 
limits) . 

In order to assess the background contribution from 
cascading, it is important to constrain as much as possible 
the template spectra and the metagalactic radiation field 
from observations. Since the cascade emission depends on 
the maximum energy and spectral index (Protheroe & 
Stanev 1993), we try to adopt realistic templates for the 
spectra obtained from observations. The total flux of the 
cascade emission also depends on the photon density of 
the MRF for which we will use models differing within 
the observational uncertainties. 

We construct a template gamma-ray spectrum for 
Eq. [3] by fitting a broken power-law model to the data 
from detected BL Lacs. For a set of j data points with flux 
dN obs /dEj(zbi) of a source at redshift zu and an optical 
depth for gamma-rays in the Universe t 77 (,2m) (10 GeV 
< E~f < 100 TeV) the flux from the same source located 
at any given redshift z = z s is 
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To test the influence of different, but realistic tem- 
plates and the metagalactic radiation field (MRF) we will 
discuss four different cases: 

1. A template obtained from observations of Mkn501 and 
1ES1959+650 and a high MRF model. These are stan- 
dard sources with similar slopes around the peak en- 
ergy. 
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2. A template obtained from observations of Mkn501 and 
1ES1959+650 and a low MRF model, to test the un- 
certainty of the MRF. 

3. A template obtained from observations of Mkn421, 
PKS2005-489, PKS2155-304 and a high MRF model. 
These sources show a steeper decline at higher ener- 
gies, which could be a sign of intrinsic absorption. 

4. A template obtained from observations of H1426+428 
and a high MRF model, to show the effect of an ex- 
tremely high value for the peak energy. 

The parameters of the HBL templates are the two spectral 
indices and the location of the maximum (the normaliza- 
tion is obtained from the luminosity function, see below). 
The numbers in Table [T] are the result of fits of the intrin- 
sic spectral energy distribution that have been calculated 
using Eq. ([8|) and the MRF model presented in Kneiske et 
al. (2002, 2004). 

For the high MRF we have assumed the "high-UV" 
model including the "warm-dust" component. The low 
MRF is the "low-IR" model without any UV emission. 
The cascade emission from an average HBL is then given 
by Eq. Q-©. The results can be seen in Fig. [2] and Fig. [3j 

In Fig. [2] the intrinsic spectra (dashed lines), the ab- 
sorbed spectra (solid lines) and the cascade emission (dot- 
dashed lines) are plotted for a source of case 1 (thin lines) 
and case 2 (thick lines). The four different panels show 
the change in the average spectrum if the HBL is located 
at different redshifts. Due to the increasing distance the 
flux is smaller and the extragalactic absorption higher. For 
comparison the data and allowed regions for the intrinsic 
spectra of Mkn501 are shown in the lower right panel at 
a redshift z = 0.03 and for 1ES1995+650 in the lower left 
panel at a redshift z = 0.047. 

In Fig. [3] the same is plotted as in Fig. [5] except for 
a source of case 3. The data and allowed regions for the 
intrinsic spectra are taken from Mkn421, PKS2005-489 
and PKS2155-304 and plotted in the panel depending on 
their redshift range. Within the uncertainties of the data 
all sources in the same figure can be fitted roughly with 
the same template. 

The location of the peak energy in the synchrotron 
spectrum of blazars apparently varies across a wide range 
from the near- infrared (LBLs) to X-rays (HBLs), although 
thermal emission from heated dust and starlight as well 
as photoelectric absorption in the EUV and soft X-rays 
hampers an unbiased measurement of the non-thermal 
spectrum. Blazars such as H1426+428 seem to exhibit, 
at least occasionally, a peak energy beyond ~ 100 keV, 
and a gamma-ray peak beyond ~ 10 TeV when account- 
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Fig. 4. The un-absorbed spectrum (dashed line), the 
absorbed spectrum (thick solid line) and the cascade 
emission (dot-dashed line) for the extreme BL Lac 
1H1428+428 at a redshift of z = 0.129. The data and 
the shaded region indicating the allowed range for the in- 
trinsic spectrum is taken from Kneiske et al. (2004). 
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Fig. 5. The intensity of the extragalactic gamma-ray 
background. The contribution from HBLs to the EGRB 
(thick solid line) for all 4 cases. (The flux without ab- 
sorption: dotted line; flux with absorption: dashed line; 
cascade contribution alone: thin solid line). Data are 
from Sreekumar et al. (1998): solid diamonds and Strong, 
Moskalenko & Reimer (2004): open circles. 



ing for extragalactic absorption (Costamante et al. 2001 
and 2003, Kneiske et al. 2004). In Fig. g] we show the ef- 
fects of cascading on such a spectral energy distribution 
(SED), adopting £ p =10 TeV, a x = 0.2 and a 2 = 1.0. The 
duty cycle of the extreme behavior is unknown, and we 
consider a sizeable contribution of such extreme BL Lacs 
as rather speculative. In Fig. Q] the same is plotted as in 
Fig.[2]except for a source of case 4 and data of H1426+428 
at a redshift of z = 0.129. 



4.2. TeV-luminosity function 

The luminosity function at gamma-ray energies (TeV-LF) 
of HBL is poorly known, since there has not been a com- 
plete survey, and the number of known sources is still 
rather low. We construct a LF based on the somewhat 
better known X-ray properties of HBLs. Bade et al. (1998) 
and Laurent-Muhleisen et al. (1999) obtained a luminos- 
ity function based on the samples from the ROSAT All 
Sky Survey, Rector et al. (2000) and Caccianiga et al. 
(2002) used samples from the Einstein Medium Sensitivity 
Survey and the Radio Emitting X-ray Sources catalog. 
For this study, we used the LF of Beckmann et al. (2003) 
who combined all the available data and derived a LF for 
BL Lacs. To obtain the TeV-LF we assume the luminos- 
ity between 0.5 keV to 2 keV L( X ) equals the gamma- 
ray luminosity above 0.3 TeV L(t c v)- This assumption 
is based on the statistics from 246 sources using an SSC 
model and relations from blazar observations (Costamante 
& Ghisellini 2002). A minimum, maximum and break lu- 
minosity of L (min) = 10 43 erg s _1 , L (max) = 10 47 erg s _1 
and L(br) = 10 45 erg s _1 respectively is used. 



5. Results 

The contribution of the HBL component for all 4 cases is 
shown in Fig. [5] Without the effect of absorption and the 
resulting cascade emission the contribution to the back- 
ground would be small. The cascade emission enhances 
the background intensity by up to an order of magnitude 
in the EGRET energy range. HBL can thus contribute 
7%, 11%, 19% or 1% to the gamma-ray background in the 
energy range 100 MeV< E 1 <20 GeV for case 1 to case 4 
respectively. 

For comparison, the contribution from EGRET blazars 
(LBL, FSRQ) is shown together with the result for case 
3 in Fig. [6l Due to the new spectral index, the change of 
cosmological parameters, the additional cascade emission 
and the reanalyzed EGRET data, the unresolved EGRET 
blazars now produce ~ 60% of the background intensity. 
Comparing the total intensity as a sum of the two con- 
tributions (thick solid line) and the EGRET data, the 
agreement is acceptable at energies below 2 GeV but is 
too small above this energy. 

Using 
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we can check the number counts of HBLs above 
300 GeV which could be observed above a flux limit of 
10~ n cm~ 2 s _1 according to our choice of the LF with 
and without extragalactic absorption (see dot-dashed and 
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dashed line in Fig. [TJ for case 1). The total population con- 
sists of about 10 6 HBLs for 0.02 < z < 5 (thin-solid line). 
We obtain 23, 24, 30 or 22 sources at redshifts z < 0.3 for 
case 1 to case 4 respectively. The peak redshift of the ob- 
served sources would be z = 0.03. The results are slightly 
higher than the number counts obtained from HBL ob- 
servations with HEGRA and Whipple. However we did 
not include the change of flux limit due to different zenith 
angles or other observational based effects. The number 
of detectable sources with the new Cherenkov telescopes 
(MAGIC, H.E.S.S., Veritas, Cangaroo) should increase by 
a factor of f observing at lower energy and in both hemi- 
spheres, which is in agreement with the ea 1 00 HBL from 
integration of the X-ray luminosity function within one 
Gpc. From Fig. 7 it is clear that independent of the tele- 
scope parameters the maximum number of sources will be 
detected at a redshift of » 0.16 (thick solid line). 

6. Discussion 

The EGRET extragalactic gamma-ray background data 
may be affected by large systematic errors. They are 
strongly dependent on the galactic foreground emission, 
which must be subtracted from the measured signal. 
Improvements in the modeling of the foreground emis- 
sion were the reason for the new determination of the ex- 
tragalactic background by Strong, Moskalenko & Reimer 
(2004). Recently, large amounts of gas, possibly connected 
with the Gould Belt, which are not accounted for in ex- 
isting HI and CO surveys, have been found using infrared 
and gamma-ray observations. This implies severe revisions 
of the gamma-ray interstellar emission models to high lat- 
itudes. Grenier, Casandijan feTerrier (2005) obtained a 
new estimate of the extragalactic gamma-ray background 
including this newly found interstellar matter. 

The largest uncertainties in determining the blazar 
contribution to the extragalactic background are their un- 
known SED and LF where only coarse estimates could 
be derived. In particular, at the low luminosity end of 
the LF, blazars could hide inside elliptical galaxies. At 
the present stage, the shape of the intrinsic spectra of 
blazars, as well as the number of blazars, seem to suf- 
fice to produce most of the background from faint, unre- 
solved sources. This has also been discussed in Giommi et 
al. (2006). They have used 5 GHz radio number counts 
and a SSC model to calculate the contribution of blazars 
to the CMB, X-ray and gamma-ray background. They 
extrapolate the number counts to very low fluxes below 
1 Jansky. As a result they found that the high energy 
gamma-ray background is overproduced by blazars if their 
SED shows an all-time high-energy gamma-ray bump. It 
is not clear if all faint radio sources which are included 
in their sample are producing gamma-rays and includ- 
ing these sources could lead to an overestimate of the 
background flux at gamma-ray energies. A more accu- 
rate treatment could start from generalizing the known 
correlations between X-ray /gamma-ray peak and lumi- 
nosity for a single, luminosity-dependent template using 
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Fig. 6. Extragalactic gamma-ray background and as- 
sumed blazar contributions due to low-peaked blazars 
(thin solid line) and high-peaked blazars (dashed line) in- 
cluding the effects of intergalactic cascading. 

leptonic and hadronic models. It is not clear whether ac- 
celerated electrons or hadrons (or both) are responsible 
for the gamma-ray emission. Both models can reproduce 
the multi-wavelength data equally well, since not enough 
simultaneous data are available (Aharonian et al. 2005). 
Therefore we have used simple power-laws for the spec- 
tral energy distribution. Recent observations at GeV and 
TeV energies have also shown flares in blazars which are 
only bright in gamma-rays. These so-called orphan flares 
have to be explained by more complex models (Reimer, 
Bottcher & Postnikov 2005). To obtain more precise blazar 
SED templates, detailed modeling of single sources with 
leptonic and hadronic models are needed, which will be 
possible if more simultaneous data from multi- wavelength 
campaigns are available. Number counts of HBLs found 
with IACTs will also be crucial to test the assumed LF 
and templates. 

Blazars show variability across a wide band of time 
scales, from minutes to years, and estimates of their con- 
tribution to the extragalactic background must take into 
account selection effects due to the excess of flaring sources 
among detected sources found in incomplete surveys, and 
on the shape of the SED. Since the duty cycle of flares 
seems to be rather small, we have ignored such selection 
effects in this work. 

Another source of uncertainty in the above calculations 
is the strength of the magnetic fields in the blazar host 
galaxies, in clusters of galaxies, and in intergalactic space. 
The magnetic field strength in the Coma cluster of galaxies 
has been measured to be (1.7±0.9) /xG (see Kronberg 1994 
for a review). The cluster magnetic field seems to be cor- 
related also with larger filaments of the matter concentra- 
tion in the Universe. Using rotation measurements of radio 
sources Ryu, Kang & Biermann (1998) could show that 
the magnetic fields in these regions could be as strong as 
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~ 1 /iG. For the gamma-ray templates considered in this 
work, the mean free path takes values between a few hun- 
dred Megaparsec and some Gigaparsecs, enough to assure 
that the gamma-rays escape freely from the intracluster 
medium before pair production occurs. In general, they 
would then enter void regions, unless the jet directions 
are correlated with the supergalactic filament orientation. 
The magnetic field in void regions is not well known, but 
it seems that it drops to 1CP 9 G or less (Lee, 01into& Sigl 
1995). Detailed cascade simulations and observations of 
gamma-ray halos can help to set limits on the extragalac- 
tic magnetic field. Aharonian et al. (2001) found no sign of 
a halo in Mkn501 in its quiescent phase. This would lead 
to much smaller magnetic fields like 10~ 16 G. The calcu- 
lations presented above were made assuming no magnetic 
field. The effects of weak magnetic fields (10~ 17 - 10- 20 G) 
on cascade emission in Mkn501 and H1426+428 are dis- 
cussed in Dai et al. (2002) and Fan, Dai & Wei (2004), 
respectively. For Mkn501 magnetic fields up to 10~ 18 G 
have almost no influence on the cascade emission at ener- 
gies above 100 McV. A magnetic field of 10 -17 G reduces 
the flux by almost one order of magnitude at 100 MeV. 
Comparing this result with the calculations of H1426+428 
the effect seems to depend strongly on the intrinsic spec- 
trum and redshift of the BL Lac. Therefore an estimate 
for the gamma-ray background seems only possible by an 
exact simulation including the extragalactic magnetic field 
as a parameter. This will be done in future work. 

Stronger magnetic fields (> 10~ 12 G) would lower the 
cascade flux from the blazars emitting beamed radiation 
by about a factor of hundred. However, the correspond- 
ingly larger number of blazar host galaxies (radio galaxies) 
would then come into play to compensate for the loss in 
flux per source. 

An analysis by Gorbunov et al. (2005) suggested a 
very weak correlation between the diffuse photons above 
10 GeV detected by EGRET and some BL Lacs (LBL and 
HBL). The low number of blazars to which the photons 
can be traced back might indicate a greater value of the in- 
tergalactic magnetic field strength, sufficient to isotropizc 
the cascade emission. 

The hard feature in the background spectrum at a few 
GeV might also indicate other emission components, with 
less broad continuum radiation. It has been suggested that 
neutralino annihilation in clumpy dark matter halos might 
naturally give rise to this bump (Elsasser and Mannheim 
2005). 

Stawarz, Kneiske & Kataoka (2005) studied an- 
other possible contribution from the parsec scale jets of 
Faranoff-Riley type I galaxies. The observed X-ray emis- 
sion in the knots of the extended jets could result in an 
inverse Compton component at gamma-ray energies. But 
neither the direct emission nor the faint cascade emission 
will lead to a larger contribution than 1%. Using a con- 
sistent model where the TeV emission is produced in the 
inner core region and in the extended jets as well could 
increase the contribution of FR-I galaxies to the extra- 
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Fig. 7. Number of BL Lacs as a function of redshift for 
case 1; thin solid line: total number of sources; thin dashed 
line: observable sources using a telescope E t hr = 300 GeV 
and a flux limit of 10~ n cm~ 2 s _1 ; thick-solid line: observ- 
able sources including the effect of extragalactic absorp- 
tion only; dot-dashed line: observable sources including 
the effect of extragalactic absorption and telescope pa- 
rameters (energy threshold and sensitivity). 

galactic gamma-ray background (for the case of M87, see 
Reimer, Protheroe & Donea 2004, Stawarz et al. 2005). 

Based on the observational finding of an increasing 
number of HBLs detected above 100 GeV energies, and 
the inescapable effect of intergalactic cascading on their 
emission, we have shown that they can indirectly con- 
tribute a sizeable fraction of the extragalactic gamma-ray 
background below 100 GeV. Ongoing observations with 
IACTs (e.g., H.E.S.S., MAGIC, VERITAS) will allow one 
to probe the assumed source population at very high en- 
ergies. Monitoring campaigns will be important to obtain 
the time-averaged flux emitted by the blazar population 
at gamma-ray energies, and to better understand the se- 
lection effects due to flaring sources. A significant increase 
in source statistics will be needed to test the assumed lu- 
minosity function. 
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